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In this paper we investigate the magnetic properties of heavy fermions in the antiferromagnetic 
and dense Kondo phases in the framework of the Kondo necklace model. We use a mean field 
renormalization group approach to obtain a temperature versus Kondo coupling (T — J) phase 
diagram for this model in qualitative agreement with Doniach's diagram, proposed on physical 
grounds. We further analyze the magnetically disordered phase using a two-sites approach. We 
calculate the correlation functions and the magnetic susceptibility that allow to identify the crossover 
between the spin-liquid and the local moment regimes, which occurs at a coherence temperature. 
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I. INTRODUCTION 

It is well known that the nature of the ground state of 
dense Kondo compounds results basically from the com- 
petition between the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) interaction and the Kondo effect. In a sim- 
ple picture it is governed by a single parameter, the 
ratio J/W, where J is the effective exchange between 
localized moments and conduction electrons and W is 
the bandwidth of the latter. The value of this ratio is 
usually tunable experimentally by pressure or composi- 
tion ratio of the compounds. The RKKY interaction is 
an indirect magnetic interaction between localized mo- 
ments, mediated by the polarized conduction electrons, 
with an energy scale of order Jrkky that pro- 

duces a long-range ordered magnetic ground state. On 
the other hand, the Kondo effect favors the formation of 
singlet states between localized moments and conduction 
electrons generating a non-magnetic ground state and, in 
the single impurity case, has a characteristic energy scale 
of order ksTx = We~^/^ . As result of the interplay 
between these two effects, some Kondo compounds are 
non-magnetic and are characterized by a heavy-fermion 
behavior (Fermi- liquid) at very low temperatures, while 
others order magnetically, generally antiferromagneti- 
cally. The study of this interplay is easily formulated 
using the Kondo lattice model, which emphasizes the 
importance of spin fluctuations neglecting charge fluctu- 
ations of the localized electrons, and has been well char- 
acterized by the "Doniach phase diagram" 0|. In this 
simple picture, the ordering temperature T^v initially in- 
creases with increasing J, then passes through a max- 
imum and vanishes at a critical coupling Jc. At this 
quantum critical point, a second order phase transition 
between an antiferromagnetic ground state for small val- 
ues of J and a dense Kondo state for strong couplings J 
occurs. This behavior of T/v has been experimentally ob- 
served in many Cerium compounds by varying the pres- 



sure applied on the system |2 Q or the relative concen- 
tration in the compounds |^;3|. 

In this paper we are interested in studying the TxJ 
phase diagram of Kondo compounds. For that purpose, 
we use an analog of the symmetric Kondo lattice with 
complete absence of charge fluctuations, the Kondo neck- 
lace model (KNM). This model was proposed by Do- 
niach 01, and its ground state has been investigated by 
a variety of methods j^-p^ . In order to assess the criti- 
cal behavior of the Kondo necklace, we apply the mean 
field renormalization group (MFRG), first proposed by 
Indekeu et al. on the KNM. This method combines 
mean field results for small clusters of spins and renor- 
malization group ideas. While mean field theory iden- 
tifies the order parameter of the cluster with the order 
parameter of the entire system, the MFRG assumes that 
the cluster order parameter rescales with cluster size . 
Using this method we obtain the phase diagram of the 
Kondo necklace as a function of temperature and Kondo 
coupling, which is qualitatively identical to the Doniach 
diagram. Within a two-site approach ||l^, we calculate 
the finite temperature magnetic susceptibility and inves- 
tigate the behavior of short-range magnetic correlations 
in the magnetically disordered phase. As a result of 
this investigation, we increase the phase diagram with 
a crossover line, which is associated to a coherence tem- 
perature that separates Kondo spin-liquid and localized 
moments regimes. 

This paper is organized as follows: In the next section 
we present the Kondo lattice and Kondo necklace models 
and discuss their analogies. In section III we apply the 
mean-field renormalization group approach in the Kondo 
necklace model and obtain the Doniach phase diagram, 
which we compare with experimental results. In section 
IV we apply a two-sites method on the non-magnetic 
phase, in order to study the role of short-range corre- 
lations on the system. We also calculate the magnetic 
susceptibility and then obtain a coherence temperature. 
In section V we summarize and discuss our results. 
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II. THE MODEL 



III. MEAN-FIELD RENORMALIZATION GROUP 



The Kondo lattice model (KLM) is a theoretical model 
for heavy fermions that can be derived from the more 
fundamental Anderson Lattice Model, in the case of well- 
developed local spin moments . It consists of two dif- 
ferent types of electrons, the localized spins whose charge 
degrees of freedom are suppressed, and the conduction 
electrons that propagate as charge carriers. It is de- 
scribed by 



H = -tJ2iclc 



h.c.) + j'^Si ■ cl,^aat3Ci,p. (1) 



The first term represents the conduction band (c| ^ is 
the creation operator, t is the hopping between nearest 
neighbors) and the second term is the interaction between 
conduction electrons and localized moments Si via the 
intra-site exchange J, where a are the Pauli matrices. 

In order to study the interplay between Kondo screen- 
ing and the RKKY interaction, Doniach proposed a sim- 
plified model related to the one-dimensional Kondo lat- 
tice, called the Kondo necklace model (KNM). In this 
model, the conduction electrons are replaced by a spin 
chain with XY coupling which eliminates charge fluctu- 
ations ||l[]: 



H 



KN 



» 3 ' 



■ Ti, 



(2) 



where and are independent sets of spin 1/2 Pauli op- 
erators. The first term mimics electron propagation, and 
in one dimension can be mapped by the Jordan- Wigner 
transformation onto a band of spinless fermions. The sec- 
ond term is the magnetic interaction between conduction 
electrons and localized spins via the coupling J, as in 
Eq. (1). 

Although the KLM is mapped onto the KNM only in 
one dimension, it is clear that, even in higher dimensions, 
Hkn has the same magnetic tendencies of the Kondo 
lattice. Since the essential features of the original model 
are kept, we expect that the main physical properties 
of the Kondo lattice will be maintained in the model 
described by Eq. (^). 

Since in the KNM approach to the Kondo lattice 
charge fluctuations are neglected, the critical properties 
in this case are described just by spin excitations. How- 
ever, it is important to emphasize, that the analysis of 
heavy fermions systems in terms of the Kondo necklace 
model is appropriate pof . A recent and very complete 
study of a heavy fermion system just at the quantum 
critical point QCP shows that a description in terms 
of local moments seems to be more appropriate for this 
kind of material. This may not be the case for all heavy 
fermions but it is also true that the universality class of 
the quantum transition of heavy fermions has not been 
determined yet and may even not be unique. 



The mean-field renormalization group was first pro- 
posed by Indekeu et al. JTst for computing critical proper- 
ties of lattice spin systems. This method has been applied 
to many statistical physics problems, both classical and 
quantum systems, with |l7|| and without disorder 
and the resulting critical exponents deviate from those 
obtained from standard mean field theories (including 
Bethe lattice calculations) JT^ . 

The main idea of the MFRG is the comparison of two 
clusters with N and N' sites respectively, subjected to 
symmetry-breaking boundary conditions. The interac- 
tions within the clusters are treated exactly, and the ef- 
fect of surrounding sites is replaced by a mean field which 
is supposed to scale in the same way as the order param- 
eter. 

For a generic spin system, one considers each boundary 
spin fixed and equal to b and b' for the N and N' spin 
clusters, respectively. After computing the order param- 
eter On and On' for both clusters one imposes a scaling 
relation between them 



ON'iK',b')^^ONiK,b), 



(3) 



with K and K' being the coupling constants of the two 
rescaled systems. Assuming a similar relation between 
the mean-fields (6' — ^6) and knowing that these fields 
have to be very small near the second order phase tran- 
sition, one can expand equation for small values of b 
and b' to obtain 



dO'j^{K',b' 



db' 



dON{K,b) 



db 



(4) 



6=0 



which is independent of the scaling factor ^. Equation (^) 
can be interpreted as a recursion relation for the coupling 
constant K, from which the critical point is extracted. 

We can apply this method to the Kondo necklace in its 
simplest version, that is, we consider two cells containing 
one and two sites each, as sketched in figure I for a 2D 
hipercubic lattice. 



a) 



b) 



FIG. 1. Clusters considered in the calculations of the 2D 
case. Each site contains two spins, one spin S (open circle) 
connected by a solid line to one spin r (filled circle). The 
interaction between sites is mediated by the spins r, where 
the dashed bonds represent interactions with the boundary 
symmetry breaking fields b' and b. 
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As we are dealing with an incipient antiferromag- 
netic ordering and the antiferromagnetism in the Kondo 
necklace occurs in the XY plane 0|, we consider a D- 
dimensional hipercubic lattice and divide the system in 
two sub-lattices A and B. The order parameter is then 
the staggered magnetization of the spins r, taken along 
the X direction (see Eq. (^) and ^ below). The x com- 
ponent of the boundary spins in the smallest cluster is 
fixed to be —h' since all the first neighbors of ri are in 
the same sub-lattice. In the two-sites cluster, ti and T2 
are in different sub-lattices, so the x component of their 
neighboring boundary spins have different signs and are 
fixed at —b and b respectively. 

Let us first consider the Hamiltonian for a one-site clus- 
ter taken on a sub-lattice A: 



Hi = J' Si • n - zb'W'rf, 



(5) 



where J' is the scaled coupling interaction, and the spin 
Ti interacts with its z nearest neighbors through the term 
zb'W'Tf. 

Similarly, the Hamiltonian for the two-sites system 
(one in each sub-lattice) is given by: 

2 

H12 = S, • n -f W{TfT^) - (Z - l)bW{Tf - t|). 



(6) 



In this case the spin ri interacts directly with 
through a term W^(Tf rf ) and both ri and T2 interact with 
their {z — \) nearest neighbors through —{z — l)bWTf and 
(z — l)bWT2, respectively. 

This method can be used to study quantum systems 
(T — 0), for which there are few available renormaliza- 
tion group techniques. For that purpose, we compute 
the ground state |0 >i and |0 >i2 of the two systems 
and their corresponding staggered magnetizations along 
the X direction. In the vicinity of the phase transition, b 
and b' can be assumed small and 

=<0|rf|0>=--^5', 
2 J 



<0|(Tf -r-)|0> 



12 



2(z-l)(Vl6j2 + 1 + 1)2 



(1 + 16 f + Vl6j2 + 1) (-1 + Vl6j2 + 1) 



where j = ^ and j' = ^■ 

The main assumption of the MFRG is the imposition 
of the same scaling relations between Mf and Mfj , b and 
b' . By doing this we arrive at the renormalization group 
recursion relation for j and j' . The associated fixed point 
equation is 



4(z-l)(Vl6j2 + 1 + 1)2 



Jc (1 + 16 j2 + Vl6j2 + 1) (-1 + ^/Wjf+T) ' 



(7) 



By solving Eq. (^ we obtain the antiferromagnetic 
quantum critical point (QCP) in any space dimension 
D. For z — 2, which corresponds to a one dimensional 
system, we obtained jc = 0.64, The critical ratios jc for 
other coordination numbers are: z — ^ ^ jc — 1-66, 
z — % ^ jc = 2.67. This values can be compared with 
recent results and the relative errors are 12.67% 

for 2D and 1.5% for 3D showing that the results are most 
relievable for higher dimensions. 

It is worth mentioning that the diagonalization of H12 
is not a simple task. However, as we are interested in 
small values of the mean field, a perturbative expansion 
can be worked out in order to obtain the eigenvalues and 
eigenvectors in powers of b [|l8| . We follow the same pre- 
scription in order to calculate the staggered magnetiza- 
tions for T ^0. 

As we have more then one variable, we cannot deter- 
mine a complete renormalization flow in the [j, T] plane. 
However, for a fixed j we can calculate T/v or vice-versa. 
In order to obtain the critical values, we consider the 
solutions of the fixed point equation associated with eq. 
(0), that is, 



9Mf(Tjv,j) 



db' 



dM(2{TN,j) 



db 



(8) 



fc=0 



where for each j < jc we obtain a TM{j) and the set of 
all critical points yields the temperature dependent phase 
diagram of the Kondo necklace. The critical line for a iD 
system [z — 6), is depicted in figure |^. Close to the zero 
temperature fixed point it behaves as: 



1^1 



\J - Jcl oc e 



N 



ksTN 



(9) 



This dependence is characteristic of this approach and 
appears for any dimension. 




FIG. 2. Phase diagram of the KNM within the mean field 
renormalization group approach (with z — 6) compared with 
experimental results for Cerium compounds extracted from 
ref. 1^^. The results are normalized by the maximum value 
of the Neel Temperature Tjv and its associated coupling J. 
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In figure II]. Endstra et al. calculated the J cou- 
pling of many Cerium and Uranium compounds based 
fundamentally on the atomic radii and interatomic dis- 
tances of these systems . More recently, Cornelius and 
Schilling, based on the variation of the lattice parameters 
generated by pressure and substitution of Si by Ge (neg- 
ative pressure) in CeM2Si2-xGex compounds (where 
M=Rh,Ru,Pd) lH-H, found a relation between the 
Neel temperature T/v and the coupling J for CeRu2Si2, 
CeRh2Si2 and CePd2Si2 H- Their results fall in 
Doniach-like curves, that when normalized by the maxi- 
mum of Tjv {T„iax) and its equivalent Jt^^^ collapse onto 
a single universal curve. 

We observe that in the experimental results, the Neel 
temperature goes to zero in the weak coupling regime at 
a value of J greater than zero, differently from the the- 
oretical predictions for the Kondo Lattice. By shifting 
the experimental results along the J axis in order to ob- 
tain Tjv = at J = 0, and using the same normalization 
used in Ref . , we obtain the results of figure ^ (2; = 6 
(31?)), where one can see that the behavior of the theo- 
retical curve is in good agreement with the experimental 
values. We can conclude that these Kondo systems be- 
have qualitatively as proposed by Doniach although in 
practice the long range magnetic order vanishes before 
the theoretical J = value is attained. 

As we are dealing with a mean-field like approach us- 
ing small clusters we do not obtain a precise value for 
the QCP for one dimension case (J = 0) such as in pre- 
vious works l^-jl^. Nevertheless, this method has the 
advantage of being very simple to implement at finite 
temperature and produce relievable values of the QCP 
for higher dimensions, as discussed before. 

The main feature of our MFRG method in its sim- 
plest approach ( 1 and 2 site cells) is that it yields the 
[T, j] phase diagram for heavy fermions as proposed by 
Doniach on physical grounds Q which, as far as we 
know, has never been obtained before. This is so be- 
cause this method captures the essential physics of the 
Kondo lattice problem, namely, the competition between 
the RKKY and Kondo interactions. However, it is im- 
portant to note that the present method can be applied 
to larger clusters, which provides better information con- 
cerning the lattice topology of the system |Q , in order 
to obtain more accurated values for the quantum critical 
point for lower dimensions. 



IV. SHORT-RANGE CORRELATIONS 

In the magnetically disordered phase, where the mean- 
fields b and b' are identically zero, we can still calculate 
some properties of the system by means of a two-site 
approximation. This method consists on solving exactly 
a Hamiltonian of two sites, i.e., Eq. (^) with & = 0. 
As the system does not present long-range (magnetic) 
order, we expect that this approach will unravel some 



aspects of this phase. We can investigate the behavior of 
short-range magnetic correlations as well as the intra-site 
correlations between localized moments and conduction 
electrons, which are related to the Kondo effect. The 
results of the calculations are shown in figures Ha) and 
1(b). 




b) 




|S|=0.1 
|S|=0.3 
|S|=0.5 



X/W 

FIG. 3. Correlation functions < S ■ r > (a) and < ri^r^ > 
(b) for different values of \5\ = \j — jc\ as a function of tem- 
perature. 

The localized-delocalized spin correlation function < 
S-T > characterizes the condensation of singlet states. At 
T = 0, this function decreases for increasing values of j, 
saturating at the value — |, which is consistent with the 
strong coupling limit where the system condenses 
into independent singlets in each site. As the tempera- 
ture increases, this correlation rapidly falls to zero and 
we observe that, since coherence in the KNM is a collec- 
tive phenomena, even low temperature excitations can 
destroy the Kondo spin-liquid regime. The inter-site cor- 
relation function < rfrf > represents the short-range 
magnetic correlations. This function has a minimum at 
a low temperature and then gradually vanishes as T in- 
creases. Differently from < S ■ t >, it has a smooth 
behavior for large T and its asymptotic behavior is inde- 
pendent of \S\ = \ j -jc\- 

The competition between these two kind of correlations 
produces a change in the behavior of the system, as illus- 
trated in figure ||. For low temperatures, the correlation 

< 5 • r > is stronger and the system is in a condensate 
of singlets, typical of the Kondo spin-liquid regime [^ . 
Near the minimum of the magnetic correlation function 

< rf rf > there is a crossover: < rf rf > begins to dom- 
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inate and destroys these singlet states giving rise to a 
paramagnetic regime with localized moments. The in- 
flection point in the magnetic correlation function, where 
it crosses < S ■ t >, defines the crossover or coherence 
temperature which separates two different regimes in the 
non-magnetic region of the phase diagram. 



0.0 



|8|=0.5 



o 
o 




<S,a:,> 



-0.8 



T/W 



FIG. 4. Competition between two kind of correlations: the 
intra-site < S ■ t > (dashed line) and the inter-site < rf rf > 
(solid line) for \S\ = 0.5. 

The magnetic susceptibility (figure |4|) , calculated ac- 
cording to eq.(p^, illustrates the difference between these 
two regimes. 



-I3E„ 



(10) 



-I3E„ 



ksT 



At low temperatures, as we are dealing with the symmet- 
ric case of the Kondo lattice, the system is an insulator 
and there is a singlet-triplet gap. At a higher tempera- 
ture there is a maximum in the susceptibility where the 
gap closes and finally, at high T, we can see an asymp- 
totic Curie- Weiss regime, typical of localized moments. 
For |5| = the system has only the Curie- Weiss regime 
and the susceptibility can be fitted as, xo « , 

The coherence temperature can be calculated by find- 



ing the zero of 
quadratic in \5\: 



and at low temperatures it is 



(11) 



The gap at T = in the critical regime, for \6\ 0, also 
has a power law behavior with the same exponent of the 
coherence temperature. 

The Doniach diagram discussed in the previous section 
can be generalized in order to include two other regimes 



in the magnetically disordered region of the phase dia- 
gram, separated by the coherence line. Figure |5| illus- 
trates the extended diagram obtained using these two 
methods. 
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FIG. 5. Complete phase diagram of the Kondo necklace 
model (for 2 = 6). The curve rjv(J) is shown here as a full 
line, whereas the dashed line represents Tcoh{J)- 



V. CONCLUSIONS 

We have successfully applied the MFRG on the Kondo 
necklace model. Using the simplest choice of clusters, we 
obtained the quantum critical point for the 2D and 3D 
cases that are in good agreenment with previous calcu- 
lations. We obtained for the first time a theoretical cal- 
culation of the Doniach diagram and compared it with 
experimental results for some Cerium compounds obtain- 
ing qualitatively good agreement. In the magnetic disor- 
dered phase we calculate, using a two-sites method, the 
short-range magnetic correlation function, the localized- 
delocalized spin correlation function as well as the mag- 
netic susceptibility as a function of temperature and char- 
acterized two well-defined regimes: at low temperatures, 
a condensate of singlets with a singlet-triplet gap and at 
high temperatures, a Curie- Weiss regime. The crossover 
line close to the QCP has been analytically calculated 
and presents a power-law behavior. 
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